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Abstract 

We present a model which supplements the Standard Electroweak Model with 
three right-handed neutrinos and one extra scalar doublet which does not develop 
a vacuum expectation value. With the aid of a discrete symmetry the neutrinos 
are kept strictly massless. This model has several interesting features. It has 
unsuppressed lepton flavour violating processes, in particular \i — > ej, hinting at 
the possibility that these may soon be within experimental reach. The Z and W 
interactions become non-diagonal at one loop level. In particular, a non-trivial 
leptonic mixing matrix is seen to arise from the clash between the charged gauge 
boson and the charged scalar interactions. 



1 Introduction 

One of the interesting features of the leptonic sector is the fact that neutrino masses are 
much smaller than the other fermion masses. To date, no experiment has unequivocally 
detected nonzero neutrino masses. In the Standard Model (SM) QXJ one precludes the 
existence of Dirac masses by not having right handed neutrino fields in the theory. 
This is a somewhat peculiar feature of the theory, since all other fermions appear as 
both right handed (RH) and left handed (LH) fields. Therefore we would like to look 
for a minimal extension of the SM which does not have this shortcoming. However, we 
must then face the problem of justifying the smallness (or inexistence) of the neutrino 
Dirac mass terms. 
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This has been done in several extensions of the SM. In some grand unified theories 
such as SO (10), a Majorana lepton number violating mass term for the right handed 
neutrinos yields neutrinos with small masses through the see-saw mechanism [0. In 
several other models neutrinos are massless and yet lepton flavour is violated. In some 
of these || ||] this is due to the lack of unitarity of the charged current mixing matrix 
which translates into non-universal couplings for the neutral current interaction. In 
others ||, R-parity violating terms in the superpotential generate flavour violating 
neutral current interactions with the electron and the u and d quarks. Most of these 
models involve rather speculative assumptions about the field content of the theory. 

On the other hand, in the SM there is no justification for the presence of only 
one Higgs doublet. In fact considerable interest has arisen in multi-Higgs doublet 
models, ever since T. D. Lee showed |J that one could have spontaneous CP violation 
in the two-Higgs doublet model. To avoid flavour changing neutral exchanges, S. L. 
Glashow and S. Weinberg , and independently E. Paschos, proposed the introduction of 
a discrete symmetry forcing each fermion to couple to only one Higgs doublet : this is 
known as natural flavour conservation ]7j . In what follows we will use a similar device 
to show how one can naturally suppress neutrino masses in a two Higgs doublet model 
with three RH neutrino fields. 

In section 2 we present our model. In section 3 we study the experimental con- 
straints arising from the processes l 2 — > l\vv, l 2 — > and Z 3 — > l 2 l\li. We also discuss 
briefly the decays Z — ► l\l 2 and W — * l\v 2 . In the last section we draw our conclusions. 
The appendix contains several derivations needed to establish the experimental bounds 
on the theory. 



2 The Model 

Our model has, in addition to the SM fields, three right handed (RH) neutrino fields [], 
N R , and one extra Higgs doublet, H 2 . The Higgs potential is chosen so that H 2 does 
not get a vacuum expectation value. In addition, a discrete symmetry is introduced 
under which 

N R -> -N R , H 2 -> -H 2 , (1) 
and all other fields remain invariant. The Yukawa and weak lagrangeans are, 

-C Y = L L H 1 T 1 C R + L L {ia 2 H* 2 )T 2 N R +h.c, (2) 
C w = ^=NL-fC L W+ + h.c, (3) 

where L L is the left handed (LH) lepton doublet, N R (C R ) is the RH neutrino (charged 
lepton) singlet and Fx, T 2 are 3x3 Yukawa coupling matrices. For later use the Higgs 
fields will be written as 

Hl = ( U v +Gr ° + ) ' H " = ( js&+ u ) ) ' (4) 



^^We could equally well construct a theory with any other number of right handed neutrinos; for 
example, one or two. We chose three since we seek a 'symmetric' theory with as many right handed 
singlet fermions as there are left handed fermions. 
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We can diagonalize the charged lepton mass matrix while keeping the charged 
current diagonal with transformations, 

Zl = (Nl,C l ) = (PlJlWI, (5) 

C R = U c Jr, (6) 

where 

Di = diag(m e , m M , m T ) = -j=uITiU Cr ■ (7) 

We can still use the fact that T 2 can be diagonalized by a bi-unitary transformation, 
T 2 = XN u Ulf R = Xdiag(ni,n 2 ,n 3 )Ul fR , and the freedom to redefine the RH neutrino 
fields by, 

Nr = Ui R u R , (8) 

to write, 

-C Y = ... + (1 + H /v)I l D 1 Ir-H~I l M^r-H + vrMI l 

+R/V2[v r Mv l + v l M^vr\ + iI/V2[v R Mv L - v L M^v R \ , (9) 

C w = -j=N Ll »C L W+ + h.c.. (10) 

Note that all fields in H 2 involve the same coupling = B^N U : the charged scalars 
link vr and II] the neutral scalars link vr and v^. 

A simple example of lepton flavour violations at low energies occurs in the decay of 
a LH lepton of flavour i into a LH lepton of flavour j and some gauge boson. This arises 
through a one loop diagram with intermediate vr and H 2 fields (Cf., for example, fig. 
1), and will be proportional to, 

3 3 

Qij = ]T M kl M) k = £ \n k \ 2 B ki B* kj , (11) 

k=l k=l 

showing a suppression in the limit in which the tik values are close to each other. 

In our model, overall lepton number conservation is imposed forbidding Majorana 
mass terms and we shall assume that there is no CP violation in the lepton sector 
(that is, the matrix N v is real and the matrix B is orthogonal; but we will keep our 
formulas general). Moreover, we shall take the new scalar masses to be within an order 
of magnitude, or so, of the electroweak scale, v = 2AQGeV. 

3 Experimental Bounds 

The first important limit on our theory comes from measurement of Gf — l/y/2v 
in the muon decay as compared to that made in the quark sector. For untagged 
neutrino flavours, we have both charged W-exchange and charged H-exchange tree 
level diagrams. Due to their different chiral structures, these do not interfere and one 
can easily find, 

Gim 5 v 4 

r(u -> euu) = ? (1 + T Vt uu Vt ee ) , (12) 

KH> ) lg27r 2 V IQMfj ' ' ^ ' 
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where M# is the mass of the charged scalars H. 

In our model, the quark mixing matrix is unitary and, therefore, the sum of the 
magnitudes squared of the elements in its first row must add to one. Any deviation 
will be a measure of how much the expression within parenthesis in Eq. ^ deviates 
from one. Using the result M, 



we find 



\V ud \ 2 + \V US \ 2 + \V ub \ 2 = 0.9977 ± .0030 (13) 
-n^n ee < 0.0848 . (14) 



v 4 



Noting from Eq. [XT] that Qu > 0, we find that, at tree level, the present sign of the 
deviation from the mean is consistent with this model. The effect of the one loop 
amplitudes will be briefly discussed at the end of this section. In any event, it is a 
good approximation to take the value of Gp as measured from muon decay, in the 
following analysis. 



The situation concerning the tau branching ratios, is still unclear |10|. This is 



particularly so for the decay into electron and two neutrinos where the relation 

T T = T»&n^fB(T^evv), (15) 
G T m T 



has only recently become consistent with the Standard Model (G M = G T ) (TT[. It is 
widely believed that this problem will continue to disappear with further experiments. 
Therefore, we shall not derive any constraints from this decay. 

As can be seen from Eq. [n], we will have lepton flavour violation in processes such 
as li —>■ Zi7 or l 3 — ► /i/^- The first of these processes has one-loop contributions from 
the diagrams in fig. 1, involving an intermediate charged scalar and RH neutrinos. 
Clearly no such diagram exists with an intermediate gauge boson. 

The most general form of the invariant amplitude for the process in fig.l is, 

M = 

= e il ui{p 1 )[im 2 <7 IJ/V k v {A 1 + £ 7 7 5 ) 

+k 2 ^{C, - D ll5 ) - + F 775 )K(p) , (16) 

where k = p — pi, gauge invariance guarantees that, 

E F 
C, = , D, = —5— , (17) 

Wl2 — mi 1712 + °m\ 

and the factor of k 2 in the 7 M terms reflects the fact that this contribution vanishes in 
the k 2 = limit. 



2 The radiative corrections are dominated by QED effects and are therefore down by a/4ir; almost 
an order of magnitude bellow our bound on new physics. However, to improve these constraints, such 
as by measuring the neutrino spectra in muon decay, these QED radiative corrections must be taken 
into account. An excellent discussion of this can be found in reference [[)]. 
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The diagrams in figs, lb and lc only contribute to the 7 M terms. So, it is easier to 
use the diagram in fig. la to find Ay, B y , E 1 and F 7 and then use Eq. [T7| to get C 7 

. We relegate the full expre' 
using Mfj » m 2 » m\, we find 



and _D 7 . We relegate the full expressions to the appendix. To lowest order in k 2 , and 



^-rar (18) 

(19) 

For the physical process I2 — > ii7, the photon is on mass shell and only the 
terms contribute, 

771 

t(i 2 ^i i1 ) = -±(\a 1 \ 2 + \b 1 \ 2 ) • (20) 

In particular, for the muon decay we can use, 

r(/i - ez/ e z/,) » « - a//) . (21) 

As pointed out above, the first sign is not a strict equality since, in this model, there 
is also an intermediate charged scalar diagram. Thus, 

B <^)^|£^« 2 > (22) 

and the model is already constrained by experiment. This should be compared with 
what would happen if only RH neutrinos and Dirac masses were added to the SM. 



In that case, to lowest order in an expansion in powers of rrik/Mw one obtains [12 
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BQ* - ery) = -£-3 | E(^) 2 £V4*e! 2 , (23) 



32tt '^.Mw 



where U is the lepton mixing matrix |TjJ and there is no term of order zero due to the 
GIM mechanism [£L4j] . Therefore, even if the neutrino masses saturate the cosmological 
bound, we get a branching ratio of the order of 10~ 40 ! Furthermore, even if one 
includes both Dirac and Majorana masses, only considerable fine tuning would lead to 
experimentally relevant values for this branching ratio ||15|| . By contrast, our model is 
already constrained by experiment. 

From the experimental result M of B(n — > ej) < 4.9 x 10 -11 we find, 



v 2 , , 



M 2 H 



Q^el < 2.8 x 10~ 6 . (24) 



It seems natural to suspect that the aren't much smaller than one, Mh is within an 
order of magnitude of v and that the B matrix might have some off-diagonal elements 
of order 0.1, say. This would raise the exciting prospect of a positive result with further 
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experiments. For the tau decays, B(t — > fi'-f) < 5.5 x 10 4 and B{r — > cy) < 2.0 x 10 4 
lead to the considerably weaker bounds, 

2 2 

-^1^1 < 22, _^|fi re |<i 3 . (25) 

Similarly, we will have flavour violation in processes such as l 2 — > 3/i due to the 
diagrams in fig. 2. Clearly, the Z diagram is very suppressed due to the Z propagator 
and we can ignore it. The general expressions have been derived in the appendix. For 
our model, 

r(/2 ^ 3/l) " 192*3 (32^)2 M^ |fi21 ' l ^ (21n( 2^ ) "l2 ) 

+i|Q 11 | 2 + l^i?e{Q 11 }] , (26) 



and, using the experimental constraints we find, 

,2 



' , fi Me | v /l + 36.4|fi ee | 2 + 2.96 J Re{fi ee } < 5.4x10" 



Ml 



^\Jl + lU\n^\ 2 + 9.28Re{^} < 90, 



2 

-^rl^el^l + 20.6|ft ee | 2 + 1.68ite{fi ee } < 48. (27) 

These are limits on different quantities than the ones showing up above, but point to 
roughly the same order of magnitude. 

In this model, one also has lepton flavour violating Z decays but these have far 
worse experimental constraints than the flavour violating lepton decays. In addition, 
the theory already predicts that, for example, B(Z — > fi + e~) should be much smaller 
than BiyH — > ey). This arises from the fact that the dominating W exchange decay 
mode of the muon has a three body final state phase space suppression, together with 
the fact that the Z has many decay modes. 

Using the results derived in the appendix we find that since Mf >> m 2 >> m 2 , 
the dominant contribution comes from the 7 M terms yielding, 

T(Z -> lj 2 ) = ^[\C z (k 2 = Mf)| 2 + \D z (k 2 = Ml)| 2 ] , (28) 

where, for M\ » Mf , 

efi?i s Ml 

From this one can easily derive a relation between the two branching ratios, namely, 

B(Z -> e-/i + ) « 3.8 x lO -5 ^ - -> e"7) , (30) 



D z (k* = Mf) « C Z (A; 2 = Mf) « ^ cot(2^) y^fr • (29) 
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showing that even for the tau (for which the coefficient is around 5.6 times larger) this 
rate is unmeasurably small. 

A similar analysis may be performed for the charged gauge interactions. The results 
can again be found in the appendix. The most important feature is the appearance 
at one loop level of off-diagonal vertices with W, leptons and neutrinos. In particular, 
the left-handed vector interaction acquires a non-diagonal leptonic mixing matrix, 

V 12 = ^f[M R ,M I ,M H ] , (31) 



where the function 



f[M R , Mj, M H \ = [~ -f— — f in (M R /M H )} + [R - /] + 1 (32) 

H Ft 

vanishes for M R = Mj = Mjj, is negative elsewhere and remains of order (— )1 for any 
reasonable ratios of scalar masses. Note that this mixing matrix is not unitary and 
therefore cannot be absorbed by a redefinition of the left-handed neutrinos. Thus we 
find the interesting feature of a non trivial leptonic mixing matrix even with massless 
neutrinos || with its possible implications for neutrino propagation in matter |4j]. 

To be exhaustive we should now go back and reassess our calculation for the muon 
decay. Indeed, the one loop corrections to the W-exchange amplitude might be com- 
parable to the H-exchange amplitude if 

^ f[M R , Mj, M h \ « J^fi^fiee (33) 

This will depend on the masses of the new scalars and on the Qy. It is also clear 
that, for example, the analysis of Beta decay must take this into account modifying 
the extraction of the CKM matrix elements. 



4 Conclusions 

We have developed a model with massless neutrinos inspired by a minimal 'democracy' 
assumption: there should exist a right-handed singlet partner for every left-handed 
particle; the fundamental scalars might also exist in several families. This was achieved 
at the expense of a 'discriminatory' discrete symmetry. 

This model has several interesting characteristics. Contrary to what happens if one 
adds massive Dirac neutrinos to the SM, in which case the cosmological limit imposes 
minute lepton flavour violations, in this model, such lepton flavour violating processes 
might be within experimental reach, especially /i — > ej. Although one also gets non- 
diagonal Z decays we found that these have levels beyond experimental verification. 
Finally, one also finds non-diagonal W interactions, and this despite the fact that the 
neutrinos are massless. A novel feature is that this comes about as the result of a clash 
between the W and the H interactions with leptons. 

This new sector of the theory will also have implications for Cosmology. In par- 
ticular, the new interactions must be weak enough to decouple the new particles early 
enough not to affect significantly primordial nucleosynthesis. This work is currently 
under way. 
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A Derivation of the invariant amplitudes 



In this appendix we derive the expressions for the l 2 — > £ l7 , l 2 — > l\Z and / 3 — > l\l\l 2 
decays, in our model. The invariant amplitude for l 2 — > £i 7 , M. = e^T^ 1 with 

= uMiim^k^ + B ll5 ) + fcV(C 7 - D, l5 ) - P(£ 7 + F ll5 )]u 2 (p) , (34) 

may be found by computing solely the Feynman diagram of fig. la. This is due to 
the fact that figs, lb and lc only contribute to the vector and axial vector parameters 
which are easier to find with the help of the Ward identity. We find 

a /7 2\ efi 2i f 1 , f xi , m 2 {xx -x 2 ) +mi(l -xi) 
A y( k = ^T~2 / dxi / dx2 XI S X2 ' 35 

OZTT 17l 2 JO JO A{Xi,X 2 ) 

) = t^t^ / dx i / rfx 2 77 v (2xi-l-x 2 ), (36) 

32n 2 Jo Jo A(xi, x 2 ) 

where, 

A(xi, x 2 ) = Mh(1 — x 2 ) — m 2 (xi — x 2 )x 2 — m\(l — X\)x 2 — k 2 [x\ — x 2 )(l — x%) . (37) 

The parameter B 1 (F 7 ) has the same expression as A 1 (-E 7 ), except for a minus sign 
for the mi term in Eq. |35] (Eq. |36[). Taking Mfj >> m\ » m\, we find the result in 
Eq. [18], to lowest order in k 2 . The parameters E 1 and F 7 do not enter in either process 
but, as mentioned above, we can use them to find C 7 and D 1 through Eq. [Tj]. The 



result, with the same approximation, is in Eq. 19 



The invariant amplitude for l 2 — > l\Z has the same structure, 

Tz = tZa(pi) [im 2 <7^(A z + £ z7s ) + 7 "(C Z - Z^ 75 ) - fc^(£7 z + F z75 )]u 2 (p) , (38) 

except that, due to the mass of the Z, the 7(U terms do not vanish in the k 2 = limit. 
The calculation is just a repetition of that for the photon. The difference in the diagram 
of fig. la (and therefore in the a^ u and k^ terms) is just due to the = cot (26 w) ratio 
between the HHZ and HH'j vertices. This is also the relevant ratio for the 7 /i 7 z, term 
of figs, lb and lc while the 7 ^ 7 _r term involves — tan 6w (and will be proportional to 
mim 2 since one chirality flip is needed in each external leg to reproduce this structure). 
It is then easy to see that a>z = ^z« 7 for a = A, B, E, F while, 

C z (k 2 ) = UvzCo + k 2 C 7 (k 2 )} , 

D z (k 2 ) = U-VzC + k 2 D^k 2 )}, (39) 



where, 



and, 



_ eQ 21 mi m 2 ri M 2 H + m 2 (t - 1) 

°° ~ 32^ mf^n 2 Jo [ M£ + m|(t - 1) J ' (40) 



rj z = tan (6 W ) tan (26 w ) + 1 . (41) 
In the same limit, M\ » m 2 » mf, we get 

eQ 2 i mim 2 



°° * 32^ 6MF ' (42) 
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A straightforward calculation leads to 
T(Z - kh) = ^[(\C z \ 2 + \D z \ 2 )-3mlRe{A z C* z +B z D* z }+mlM 2 z /2(\A z \ 2 +\B z \ 2 )} , 

LZ7Y 

(43) 

with the parameters evaluated at k 2 = M§. For M\ » M| we can use the ap- 
proximate expressions derived for the photon parameters in lowest order of k 2 /Ml to 
get, 

C z (k -M z ) - -—^ z{ -— + -— + (M z /M H )} 
ett 2 i M 2 Z 

and similarly for D z . 

A similar calculation shows that the W interactions also become off diagonal at one 
loop level. For the Z and 7 interactions, the off-diagonal couplings are excluded by 
symmetry at tree level. Therefore, they are protected against one loop divergences. For 
the W interactions the situation is different since the lack of such off-diagonal couplings 
corresponds to a basis choice and is not dictated by any symmetry. In this case, the 
infinities are cancelled by the contributions from the counterterms, which we calculated 



in the on-mass-shell renormalization scheme JL6). For the process 1 2l — > W v\ h , we 
find, 



w 



u Vl { Pl )[im 2 a^k u A w {l + 75) + YC W {1 - 75) - WE W {1 + 7s)Mp) , (45) 



where, in the limit that the square of the new scalar masses (Ml for the charged 
scalars, M\ for the neutral scalar and M 2 for the pseudoscalar) are much larger than 
M^, and to first order in k 2 , we get, 

Aw = ^= |^{[ 6(M 2 l _ M 2 h? (M 2 H ~ M 2 R + M 2 R In (M 2 R /M 2 H ))} + [R -> /]} ,(46) 
C « = + - il + 1> . (^) 

Ew = ^ 3^ m2(1 3(M^M^)3 ( 2(M ^ ~ Ml) + + Ml) ln ( M ^0 )1 

(48) 

As an easy check, we note that, except for a factor of g/\/2 instead of e, we recover 
the results for l 2 — > £17 when Ml = M 2 = Mfj. In particular, in that limit, we obtain 
Cw = reflecting the fact that the vector coefficients for l 2 — > ^7 vanish at zero 
momentum transfer. 

For the process l 2 — > hhh, we get contributions from fig. 2. Note that there is 
no WW-box contribution since, at tree level, there is no mixing matrix in the lepton 
sector. Similarly, chirality considerations exclude the WH-box diagram. 

The amplitude for this process can be written as, 

M{1 2 -> hhh) = M( Pl , P2 ) - M( P2 , Pl ) (49) 
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with, 

M( Pl , P2 ) = M^( Pl , P2 ) + M H ( Pl , P2 ) + M z ( Pl , P2 ) (50) 

where M. 1 is the one photon exchange amplitude and Ai H is the HH-box contribution. 
J\A Z is the Z-exchange amplitude which is very suppressed by the Z propagator since 
the momentum transfer is bound by m|. Therefore we will ignore this contribution. 
Using the expression found above, it is easy to get, 

•M 7 (pi,p 2 ) = [T^]Mp 2 )e 7 ^3)]Ai , (51) 
with ki — p — Pi- Similarly, 

M H ( Pl , P2 ) = S'[«i(pi)y t (l - 75)« 2 (p)][«i(P2)7m(1 - 75)ui(p 3 )] , (52) 
where, in our model, 

In the limit m 2 >> m 1; this leads to a decay rate, 

r(/2 ^ 3/l) = 768^ (r7 + ^ + r7H) ' (54) 

where 

= 4[4ln{^)-^-]{\A\ 2 + \B\ 2 )-12Re{AC* + BD*} + 3{\C\ 2 + \D\ 2 ) , 

= iq\s\ 2 , 

n H = 8Re{{C + D-2A-2B)S*)} , (55) 

and the subscript 7 is implied. 

A process such as Z 3 — > Z1Z1Z2 has a similar expression in the limit m 3 >> m 2 , m b the 
only difference showing up in the fl factors and in the logarithm. The processes with 
no identical particles in the final state are, of course, easier to calculate. They lead to 
limits on different combinations of fl factors but give no qualitatively new information. 
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by the United States Department of Energy, under the contract # DE-FG02-91ER- 
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FIGURE CAPTIONS 



Figure 1: One-loop contributions to the lepton flavour violating l 2 — > Z17 vertex. 
The diagrams in figures lb and lc only contribute to the vector and axial-vector pa- 
rameters. 

Figure 2: One-loop contributions to the process I2 — > 3Zi. In our model, there are 
contributions from a penguin diagram (fig. 2a) and from a HH-box diagram (fig. 2b). 
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